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Abstract 

A general formula for the current through a disordered normal- 
superconducting junction is derived, which is valid at finite temperature and 
includes the full voltage dependence. The result depends on a multichannel 
scattering matrix, which describes elastic scattering in the normal region, and 
accounts for the Andreev scattering at the NS interface. The symmetry of 
the current with respect to sign reversal in the subgap regime is discussed. 
The Andreev approximation is used to derive a spectral conductance formula, 
which applies to voltages both below and above the gap. In a case study the 
spectral conductance formula is applied to the problem of an NINIS double 
barrier junction. 
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I. INTRODUCTION 



The study of electronic transport in normal-metal - superconductor (NS) or semiconduc- 
tor - superconductor (SmS) sandwiches has attracted a considerable amount of interest in 
the past years. At sufficiently low temperature and in high quality mesoscopic samples, the 
phase-breaking length of the electrons is larger than the typical system size. The resulting 
coherence of the electron wavefunctions produces directly observable quantum effects, which 
manifest themselves in a typical non-local response of the currents to the applied fields. Of 
special interest is the effect of the electronic phase-coherence in a normal metal - supercon- 
ductor system. In the standard theory of the proximity effect, the influence of the supercon- 
ductor on the normal metal can be understood in terms of the coherent coupling of electrons 
and holes in the metal as described by the Bogoliubov-de Gennes (BdG) equationsS. The 
correlation between electrons and holes is produced by the process of Andreev reflection!'! 
at the NS interface, which mixes electron and hole states while quasiparticle current is con- 
verted to supercurrent. This microscopic picture of electron-hole correlation is equivalent 
to that provided by a condensate wavefunction in the normal metal which is induced by 
the superconductor through the continuity conditions at the NS interface. The scattering 
matrix approach makes use of the microscopic single particle picture of coupled electron 
and hole channels providing a straightforward and powerful tool within a formalism of the 
Landauer-Btittiker typJHl. 

The quality of the interface as well as the phase breaking processes determine the strength 
of the proximity effect and naturally have their impact on the current-voltage character- 
istics (CVC). A few fascinating transport experimentsIHl have been carried out recently, 
investigating temperature and voltage dependence, as well as the flux modulation, of both 
NS and SmS junctions. Interestingly, the relative strength of the interface barrier and the 
elastic scattering in the normal region is crucial for the features of the CVC. The ratio of 
the two determines whether subgap conductance peaks arise at zero or finite voltage0Hli. 
The investigation of these so-called zero and finite bias anomalies in the subgap conduc- 
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tance have been the object of recent experimentsEjila. The present work draws much of its 
motivation from the ongoing discussions and experiments in this area. In our case study of 
a double barrier NINIS junction we observe zero and finite bias anomalies and shed light on 
the mechanism producing these structures. 

Stoof and NazarovEl have recently described carrier transport in disordered NS junctions 
in terms of an energy- dependent diffusion constant, successfully explaining the recent ex- 
periments on reentrance in the conductivity at low temperatures!'!. Their work uses the 
quasi-classical Green's functions technique, which allows to describe transport both close to 
equilibrium^ and far away from equilibrium0@ and facilitates the averaging over disorder 
in diffusive conductors. This approach is quite indispensable if phase-breaking processes are 
to be included. An appealing alternative approach is the scattering matrix technique, which 
relies on the quasiparticle wavefunctions described by the BdG equations. While being valid 
in a general context, it describes the transport in mesocopic systems from a ballistic point 
of view, suggesting an intuitive, nearly classical pictureBi. The transport through normal 
or superconducting leads is expressed through the properties of a multichannel scattering 
matrix accounting for all elastic scattering processes, whether they be due to a geometric 
constriction, single impurities, or disorder (inelastic processes are excluded from such a de- 
scription). By these means, the transport problem is reduced to solving a ballistic problem 
at the interfaces of the normal and superconducting leads. The current and the conductance 
of the system are determined analytically in terms of the transmission and reflection ampli- 
tudes of the scattering matrix. Adhering to this formalism rather than the Green's function 
technique helps us to improve our understanding of zero and finite bias anomalies. 

The study of normal-superconducting junctions goes back to the works of Kulik@ on SNS 
junctions and of Blonder et al.i on NIS junctions, who studied nonlinear transport within the 
framework of the BdG equations considering quasi one-dimensional models. On the other 
hand, the scattering matrix technique was developed by Landauer and co-workersi@ in the 
linear response regime, resulting in the well known conductance formula for a normal metal. 
Lambert^ and Takane and Ebisawail extended the approach to include superconducting 



segments, on the basis of which Beenakkera derived a zero temperature, linear response 
conductance formula for the transport through NS junctions. A few studies have been 
carried out recently0@ combining the scattering matrix approach with the finite voltage 
transport model of Reft They have limited themselves to the energy dependence of the 
scattering states to extract finite voltage properties of the CVC in the subgap regime. Here, 
we extend these works to voltages above the gap and additionally take into account the full 
voltage dependence of the transport problem. This is important within the context of the 
sign reversal symmetry of the differential conductance which we discuss below. At the same 
time, we provide a common framework for the above studies, tracing them back to a single 
general formula. 

In the present paper, we derive the general expression for the current through a NS 
junction in the scattering matrix approach, valid for multiple channels, finite voltage, and 
nonzero temperature. In section II, we present the derivation of the current- volt age relation 
and express it in terms of a spectral conductance formula, thereby accounting for the full 
voltage dependence of the transport problem. We discuss the related symmetry of the CVC 
with respect to the sign reversal of the bias, which follows independently of any approxima- 
tions for the BdG equations. In section III, the reflection at the NS interface is made explicit 
using the Andreev approximation and a spectral conductance formula is obtained express- 
ing the result in terms of the normal, energy and voltage dependent scattering matrix. We 
illustrate this formula in the new regime of voltages above the gap and expose its connec- 
tion to previously obtained limits. In section IV, the conductance formula is used for the 
analytical as well as the numerical analysis of a multichannel, double barrier NS junction. 
We describe the existence of resonances due to quasi-bound Andreev states, and show that 
they produce sharp conductance peaks in both the single and multichannel junction. We 
explain the generic mechanism underlying the appearance of zero and finite bias anomalies 
in the ballistic two barrier system. Furthermore, we interpret our results in connection to 
experiments and theories of zero and finite bias anomalies in disordered NS junctions. 
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II. CURRENT VOLTAGE RELATION AND SPECTRAL CONDUCTANCE 



We consider a normal-superconducting junction with quasi one- dimensional, ballistic 
normal and superconducting leads, as shown in Figure 1. The pair potential vanishes in the 
normal part, due to the absence of attractive electron-electron interactions. Between the 
normal lead and the NS interface, the electrons traverse a disordered region, the transmission 
through and the reflection thereof are described by a scattering matrix. The disorder stands 
for any source of elastic scattering processes, including effects of a magnetic field. Between 
the scattering region and the NS interface, a small ballistic normal region serves to separate 
the scattering in the normal part, which mixes all electron channels at a given energy, 
from the scattering at the NS interface, where electron and hole channels are mixed in the 
reflection process (the evanescent modes in this intermediate ballistic region are neglected). 
The coherent scattering in the disorder region of the normal metal is described by the 
AN x AN scattering matrix, 



which we denote by Sn- The matrix connects the incoming N electron (hole) channels 
If (if) on each side to the equal energy outgoing channels 0\ (Of) according to Figure 1 
(i = 1,2). The N channels represent the different transverse states at the Fermi surface 
(we neglect a change of the number of channels with increasing voltage, which is strictly 
valid only in the presence of electron-hole symmetry^). r# and are N x N reflection 
and transmission matrices for electron channels, fa and Uj the comprehensive 2N x 2N 
matrices including the complex conjugated reflection and transmission amplitudes for holes. 
(Here we will denote the complex conjugate of a matrix m by m*, the transposed matrix by 
m T and the adjoint by rrv .) Following usual convention, we include the propagation in the 
ballistic region N 2 in the scattering matrix. For states normalized to carry unit probability 




(1) 
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currents, the continuity equation implies the unitarity of the scattering matrix. Note that 
the important distinction to previous work is, that we allow the scattering matrix not only 
to depend on the energy of the states, but also intrinsically on the applied voltage (the 
voltage dependence is not explicit in the notation). This is necessary to account for the full 
voltage dependence of the scattering problem. We can thus describe the deformation of the 
states due to the space dependent potential including, e.g., the voltage dependent Schottky 
barrier at a semiconductor-superconductor (SmS) interface. 

We define an analogous unitary scattering matrix for the NS interface Si by 
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Incident and outgoing channels are again labeled according to Figure 1. r, r' and t, t' are 
N x N reflection and transmission matrices between states normalized to unit probability 
current in the normal and superconducting lead0, and are grouped into the 2Nx2N matrices 
f/, r'j, tj, t'j. In the evaluation of the current, we will consider the matrix ([]]) to be specified 
by an arbitrary model for the disorder, while the matrix (|2]) will be made explicit below 
using the BdG equations. For the present purpose they are both arbitrary and thus the 
shape of the pair potential need not be specified. The total effect of all scattering processes 
in the disorder region and at the NS interface can be described by a global scattering matrix 
of the type (|2|), which is unitary as well. We restrict ourselves to its submatrix R describing 
the reflection into the normal region, 
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R ee , R e h, Rhe, and Rhh are again N x N reflection matrices. R (e, V) will be computed below 



from the given scattering matrices ([TJ) and (0). We now derive the current-voltage relation 
based on this global reflection matrix (we do not need the transmission matrix, since we 



determine the currents in the normal lead). 

Applying a voltage V (denoting the voltage in a two point measurement) on the normal 
side has two consequences. First, the voltage induces an electrostatic potential drop over 
the disorder region in the NS junction, resulting in a deformation of scattering states. The 
coupling of incident and outgoing channels is thus voltage dependent in general, as described 
by R(e, V). A stationary state incident from the normal lead (of energy e in channel v) 
consists of the incident electron and the reflected electron- and hole- states and carries the 
current (e = |e|), 

h (e, V) = =^ (l-E I Ree (6, V)« v | 2 + £ I RHe V) „ v \*\ . (4) 
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Second, the applied voltage shifts the chemical potential of the reservoir attached to the 
normal lead by — eV with respect to the reservoir on the superconducting side. The de- 
formation of the states by itself produces no net current!!. The net current flow results 
exclusively from the difference in occupation of the (finite voltage) scattering states incident 
from the left and right reservoirs. Writing the sum over channels as a trace, we obtain the 
current- volt age relation 

/ = Jde±[f(e)-f(e + eV)]G a (e,V), (5) 

with the spectral conductance 

G s (e, V) = ^Trjl - R\ e (e, V) R ee (e, V) + R{ e (e, V) R he (e, V) }. (6) 

A factor two accounts for the spin degeneracy of the channels. The defined spectral conduc- 
tance G s (e, V) describes the current contribution of the incident scattering states at energy 
e, at a given voltage V (by convention, the energy is measured with respect to the chemical 
potential in the superconductor). Formulas (|3p and (((]) imply the differential conductance 

dl 



n r 

W \v = -] def'(e + eV)G,(e,V) 



+ tdi i[f ie )-J {e + eV)]^hn, (7) 



with the expansion 

J£ \ v = G s {-eV, 0) + 2Vd v G s (e, 7) | £ =- e v,v=o +••• , (8) 

at zero temperature. This differs from the differential conductance of Ref.i, dl/dV = 
G s (— eV, 0), by accounting for the change in the conductance of the open channels with 
increasing voltage. 

To complete the general derivation we need to express the matrix R in terms of the given 
scattering matrices (Jl]) and (0). Summing over all scattering paths of an incident electron 
or hole excitation, multiply scattered between the disorder region and the NS interface, we 
arrive at 

R (e, V) = hi (e) + k 2 (e) [l - 77 (e) f 22 (e) ] "V, (e) hi (e) . (9) 

Apart form the direct reflection at the disorder region, the simplest process contributing 
consists of an excitation, which is first transmitted (£21) through the disorder region, reflected 
(77) at the NS interface, and transmitted (ti 2 ) back to the normal lead. All further paths 
result form iterative scattering processes between the disorder region and the NS interface 
(note that the scattering matrix for the disorder includes the propagation in the ballistic 
lead N 2 ). The expressions (|^), (|B|), and ([]) determine the general form of the current- 
voltage relation of a disordered NS junction, without having made any assumptions about 
the nature of the scattering at the NS interface like, e. g., the shape of the pair potential A. 
In the next section, we evaluate the spectral conductance (§) further by using the Andreev 
approximation for the scattering at the NS interface. 

We close this section with a discussion of the symmetry of the CVC with respect to sign 
reversal of the applied voltage. In the subgap regime e\V\ < A, the incoming quasipar- 
ticle excitations may not enter the superconductor. The probability current of the states 
|e| < A is totally reflected and thus the global reflection matrix R(e,V) of (||) is uni- 
tary. The unitarity produces the relations R\ e R ee + R/ he Rh e = 1 and R ee R\ e + R e hR\ h = 1. 
The symmetry of electron- and hole- type excitations in the BdG equations guarantees 
Reh (e, V) = —R* he (— e, V). As a consequence, the subgap conductance takes the form 
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G s (e, V) = ^"Tr{4 e (e, V) R he (e, V) } = ^{4h V) ^ (e, V) } 

= ^Tr|4 e (-e, V) (-e, V) } = G s (-e, V) . (10) 

A subtle issue is that this symmetry does not yet imply a symmetry in the CVC under 
reversal of voltagell. The latter requires that G s (e, V) = G s (— e, —V), which amounts to 
G s (e, V) being independent of voltage. Then we have G s (e) \ e =- e v = dI/dV\y and the 
differential conductance is invariant under sign reversal of the voltage. Indeed, in recent ex- 
periments on SmS junctionsEzE, an asymmetry in the CVC was found in the subgap regime, 
which can be understood on the basis of the above discussion taking into account the voltage 
dependent Schottky barrier at the SmS interface. The deviations from the symmetry are of 
the order of eV7/i or eV/V , where \i denotes the chemical potential and V characterizes the 
strength of the scattering potential. An explicit account of the voltage dependence of G s 
requires the scattering matrix Sn to be determined in the applied electrostatic potential. In 
principle, this task demands the self-consistent solution of the scattering problem and the 
Poisson equation!^, a problem which has not yet been thoroughly addressed. In several cases 
of interest, though, an approximate consideration of the voltage dependence will furnish an 
accurate description. 



III. SPECTRAL CONDUCTANCE IN THE ANDREEV APPROXIMATION 

We start from the spectral conductance formula @ and evaluate it by solving the bound- 
ary conditions for the transparent interface in the Andreev approximation. The stationary 
states in the ballistic leads are solutions of the BdG equationsB and are of the plane wave 
type. A step function model for the pair potential A (x) = A Q e lx 6 (x) is assumed, which 
neglects the suppression of the pair potential in the superconductor on the distance of a 
coherence length. The NS interface connects electrons and holes of the same channel with 
a reflection amplitude depending on the reduced chemical potential \i v = /i — h 2 k 2 L /2m. In 
the limit e, A «C fi u , the BdG equations are simplified by linearizing the dispersion relation 
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around the effective Fermi wave number = ^2m^ v /h. The boundary conditions are 
fulfilled in the Andreev approximation, which treats the wave number k v in the phase factor 
of the single-particle excitations only to zeroth order. As a consequence, incoming electrons 
are purely reflected into holes and vice versa. The reflection of electron (hole) channels into 
hole (electron) channels is thus described by scalar reflection amplitudes, which turn out to 
depend only on the energy e and not on the channel index. The reflection matrix at the NS 
interface is given by 

^ n (A \ ( n „-«xr (A \ 



77(e) 
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The 2N x 2N global reflection matrices R ee and Rh e can be determined from equation 
and using (Bh, we obtain the multichannel spectral conductance formula 
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G s (e,V) =?L(l+|r(e) | 2 )Tr{4(e) 

x (l- I r (6) I 2 r 2 T 2 (- e ) r^ 2 (-6) ) [l - T (e) 2 r 22 ( e ) r 2 * 2 (- e )] _1 1 21 ( e )} , (13) 

valid at all energies. Eq. (|13|) is one of the central results of this work. Combined with 
equation (||), it provides the finite voltage, finite temperature CVC of a disordered normal- 
superconducting junction in the Andreev approximation. The spectral conductance depends 
on the scattering matrices of the electrons at energies ±e as a signature of the presence of 
Andreev reflection (the explicit voltage dependence, although present in both reflection and 
transmission matrices, is not indicated in the notation). The dependence of this formula on 
the phases of the reflection and transmission amplitudes proves crucial in determining the 
resonance peaks in the conductance. The elementary process, which contributes to these 
phases is the propagation of an electron and a hole between the disorder region and the NS 
interface. 
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If no inter-channel mixing takes place, i. e., the matrices tij and ru are diagonal, the 
conductance reduces to the quasi one- dimensional form, 
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TV = ir 22) vv is the reflection amplitude for a state coming in from the right side of the 
scattering region (Fig. 1) and R u = \r u \ 2 and T v = 1 — R v denote the reflection and 
transmission probabilities of the Z/'-th channel. The last term in the denominator describes 
the crucial scattering process which involves twice the propagation between the disorder 
region and the NS interface, once as an electron and a second time as a hole. 

For voltages well above the gap, \e\,V ^> A (still assuming |e| /i), the Andreev 
reflection is strongly suppressed and drops according to V (e, V) ~ A/2|e| — >• 0. The spectral 
conductance (0) asymptotically approaches the expression for a normal junction, 

G s (6, V) = ^Tr {4 (e, V) t 21 (e, V)} (15) 

(note that this Landauer type formula is the spectral conductance at finite voltage and 
energy and deviates from the differential conductance). The conductance ( |15D exhibits no 
particular symmetry properties with respect to sign reversal of the bias. 

At voltages below the gap, we make use of | T (e) |= 1 and T v = 1 — R v and obtain the 
spectral conductance (0) in the form, 

W ) = ^i; TA,,v)TA r ,v) 

h Sil+fi, (e, V) R„ (-6, V) - 2Re [V (tf r„ (e, V) r"„ (-6, V) 
The reflection and transmission coefficients at ±e are symmetrically involved in this for- 
mula, which results in the symmetry of the CVC discussed in section II. In contrast, the 
spectral conductance (0) at voltages above the gap becomes increasingly asymmetric as it 
asymptotically approaches the Landauer expression. As an important difference to the nor- 
mal expression ([15]), the conductance ( |T6| ) contains the phase information of the scattering 
processes, which is imperative for the distinction of zero and finite bias peaks in the double 
barrier NS junction dicussed below. 
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The linear response limit (e, V — * 0) of ([16]) can be determined using T (oy = 1 and 
R v = 1 — T v and takes the form! 

q (q\ _ ^ T u (0) , . 

1 ; h ^(2-T,(0)) 2 ' 1 ; 

Eq. (|17j) has proved very useful through its remarkable simplicity. 



IV. DOUBLE BARRIER NS JUNCTION 

In this section we apply the above results to a double barrier NINIS junction, which is 
a model system for the study of the interplay between normal- and Andreev- levels in a 
Fabry-Perot type IiNI 2 interferometer. We will learn about the mechanism govering the 
appearance of zero and finite bias anomalies and give insight into the interpretation of the 
corresponding anomalies in a dirty NS junction. First, we discuss the structure in the 
conductance of a single channel NI1NI2S junction, which we trace back to the presence of 
Andreev resonances. Secondly, we present numerical results for a multichannel junction, 
showing that the typical resonance structure of a single channel survives the summation 
over the channels. This stability is a peculiarity of the superconducting system not observed 
in a normal NI1NI2N double barrier junction. 

Since the channels separate in the double barrier problem, we can make use of the result 
( [TID for the conductance gE. G s depends on the phases (p (±e) of the reflection amplitudes 
r (±e) as well as on the complex amplitude T (e) of the Andreev reflection. We use the 
notation T (±e) = T±, R (±e) = R±, and r (±e) = y/R±e l(p ( ±e ' for the reflection amplitude, 
the phase factors being determined by the potential barriers I\ and I2 and the propagation 
between them (for simplicity, we omit the explicit dependence of the scattering matrix on 
voltage). We rewrite the amplitude of the Andreev reflection as Y (e) = |r|e _ "^ e ) with the 
phase $ (e) = arccos (e/A) below the gap and vanishing above. The conductance simplifies 
to 



2(1 + |r| 2 )e 2 








T+ (1 - 


r| 
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(18) 
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which is always less or equal to the universal value 4e 2 /h. Note that the Andreev reflection 
is supressed above the gap (|T| < 1), while the phase •& (e) vanishes. 

Below the gap, the phase $ (e) is decisive for the resonances produced by the interference 
of multiple scattering processes. The maximal conductance (which is twice the normal 
conductance 2e 2 /h) is only assumed if the reflection probabilities R± are equal and the 
phase if (e) fulfills the resonance condition, 



The analogous property is found in a double barrier system within a normal constriction 
(NI1N2N), where the maximal conductance 2e 2 /h is reached if the reflection probabilities of 
the two barriers coincide at the resonance energy. 

Using the conductance formula ([18|), we consider a one channel NINS junction, which 
consists of a ballistic NS junction containing a single potential barrier at a distance d from the 
perfect NS interface. In the limit of a high potential barrier, R + and i?_ are approximately 
equal. The reflection amplitudes r (±e) describing the propagation of electrons and holes 
have nearly constant modulus and phases ip (±e) = it + 2k±d. Using the wave number 
k± = mvp ± e/vp (vp is the Fermi velocity of the channel), the resonance condition ( |T9"D 
yields the spectrum of Andreev levels, 



which predicts resonances in the conductance of a typical width proportional to the trans- 
mission T of the barrier (similar Andreev resonances are found for voltages above the gap 
at e n = nirvp/2d with a width roughly oc T/T(e)). The phase •& (e) varies between —tt/2 
and from e = to e = A and guarantees the existence of at least one Andreev resonance 
for arbitrarily small d. In the d — > limit, the weight of this resonance lines up with the 
gap voltage and we recover the NIS junction as discussed by Blonder et al.i. featuring a 
suppressed subgap conductivity and a peak in the differential conductance at the gap volt- 
age. This peak can be understood in terms of the Andreev resonance which moves to the 
gap energy for d — > 0. 



cos [<p (e) - (p (-e) - 2$ (e)] = 1. 



(19) 




(20) 
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We now introduce an additional barrier at the NS interface and analyze the resulting 
double barrier NI^N^S junction using the same conductance formula flI8|) . Following our 
definitions, <p (±e) represent the phases for the reflection of electrons entering the double 
barrier scattering region from the right (superconducting side). The corresponding reflection 
amplitudes are given by 

r(±e) = r2 + T^r^' ( 21 ) 

where ri, U are the amplitudes of the left (i = 1) and right (i = 2) barrier. The phase of this 
reflection amplitude plays the major role in determining the structure of the conductance, 
as it controls the existence of resonances according to (fL9|) . Let us fix the barrier Ii and 
increase I2 slowly, keeping their strengths I x > I 2 . In this situation, the INI interferometer 
develops pronounced Andreev resonances. For r\ ^> r 2 , the phase ip (±e) of the reflection 
amplitude r (±e) ~ t2r\e 2lk±d produces a linear energy dependence of the phase ip (e), which 
changes by 2n on the scale vp/d and results in equidistant resonances^, in accordance with 
(p0|). As the strength of I2 is increased, the resonances pair up as is illustrated in Figure 2. 
The phase function ip (e), as displayed in Figure 2 (solid line), can be used to determine the 
location of the resonances by finding those combinations of energies ±e, which have a phase 
differenceii Atp (e) = (p(e) — ip (— e) = it + 2nn. The double period of A(p (e) with respect 
to the period of ip (e) accounts for the pairing of the resonances. 

When the strengths of the barriers become of the same order, Ii ~ I2, the spectral weight 
of the INI interferometer is shared by Andreev quasi-bound states of a mixed electron-hole 
character and normal electron quasi-bound states. Due to the large gradient of the phase 
close to the normal resonances, see Fig. 2, the Andreev resonances tend to be pinned to 
normal resonances at either +e or — e (a notable exception to this rule is found when the 
normal resonance is aligned with the Fermi energy, but the Andreev resonances remain at 
finite bias). While the Andreev bound states contribute to the current, normal bound states 
do not couple to the superconductor and thus do not participate in the charge transport. 
This is reflected by the symmetry of the G s (e) under reversal of voltage which is is observed 
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for all barrier strengths, see Figure 2. 

As the barrier strength is increased further to I2 > Ij, the Andreev resonances are 
weakened and eventually disappear. Although the normal resonances dominate in the INI 
interferometer in this regime, only the weak Andreev resonances show up in the conductance, 
which thus exhibits only a weak (but still symmetric) subgap structure. The phase function 
ip (e) in Figure 2 (dashed line) becomes nearly constant for r\ <C r 2 , see (pl~l), and the phase 
condition for resonance (|19[) cannot be met. 

Let us compare the transport in the double barrier systems NI1NI2S and NI2NI1S, i.e., 
with inverse sequences of the barriers Ii and I 2 . Note that the transmission T (e) is identical 
for both cases, see Figure 2 (dotted line). Both the (finite voltage) conductance in the normal 
NININ junction ([TJ) as well as the linear response conductance in the superconducting NINIS 
junction (|T7| ) are thus independent of the sequence of the barriers Ii and I 2 . Let us assume 
that Ii ^> I2. In this first barrier sequence, we have a strong energy dependence of the phase 
(p (e) (solid line in Fig. 2), which implies the existence of Andreev type resonances at finite 
bias. The electrons entering the INI interferometer from the normal lead, are given enough 
time to build up an Andreev resonance and preferably leave into the superconductor. For 
the inverse barrier sequence, the barrier I 2 at the NS interface dominates. The weak energy 
dependence of the phase ip (e) of reflection allows no sharp resonances to build up (dashed 
line in Fig. 2). This reflects the fact that the electrons which enter the INI region leave 
through L into the normal lead before an Andreev resonance can build up. In summary, the 
spectral density in the INI interferometer changes radically with the coupling strengths of the 
normal and superconducting leads. Normal resonances dominate when the interferometer is 
coupled more strongly to the normal lead, whereas Andreev resonances take over in weight 
when the coupling to the superconductor is stronger. At any instant, however, only the 
Andreev states participate in the charge transport. 

We turn to the numerical analysis of a multichannel NINIS junction which we carry out 
using the conductance formula (13). This formula allows to extend the linear response study 



of Ref.E3 to finite voltage and temperature. We investigate an NI1NI2S junction with two 
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^-function barriers of typical strength H = J V (x) dx ~ Tivf and corresponding reflection 
probability R = H 2 / (^H 2 + h 2 v F ^j assuming values between R = 0.2 and R — 1. We vary 
the relative barrier strengths to cover the range between the two limits Ii > I2 and Ii < I2 
discussed above. The distance between the barriers is chosen to be of the order of or larger 
than the coherence length of the superconductor, providing the forward channel with one 
to a few Andreev resonances. The number of resonances increases with the incidence angle 
of the channels. We chose leads with a cross section area of (100/ kp) 2 , which amounts to 
about 800 transverse channels. The ratio of the energy gap to the Fermi energy is assumed 
to be A/eF = 0.002. Each channel features the typical conductance structure of paired 
Andreev resonances exposed above. Their positions and widths depend on the ratio of the 
barrier strengths Ii and I2 as well as the longitudinal kinetic energy of the single channels. 
Remarkably, the overall conductance, which is obtained through summation of single channel 
conductances, exhibits a characteristic subgap structure signalling the presence of Andreev 
resonances. In contrast, the overall conductance of the corresponding NI1NI2N junction is 
practically constant, the normal resonances of the INI region having been averaged out. 

The numerical study of a 3-dimensional NINS junction has shown that both the positions 
and the number of resonances in the overall conductance correspond to those in the forward 
channeflii. In the NINIS junction, we do not find a direct correspondence of the resonances 
of the total conductance with the forward channel nor with any other specific channel, 
although a clear resonance structure still survives the summation over the channels. 

Let us concentrate on the conductance formula (|I"6D, valid at subgap voltages, and on the 
properties of the CVC close to zero voltage. For I x > I 2 , the denominator of ([16]) changes 
rapidly with the strong energy dependence of the phase of r„ (e) which is responsible for 
the appearance of conductance peaks at finite voltage. The pronounced structure in the 
conductance survives the summation over the channels as displayed in Figures 3 and 4 
(solid lines). The repulsion of the Andreev levels around zero voltage produces a minimum in 
dl/dV at zero voltage. For Ii < I 2 , the phase of the reflection amplitudes r v (e) has negligible 
energy dependence and the numerator of (pjl) dictates the features of the conductance. The 
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expansion of the product T u (e) T v (— e) = T 2 — e 2 T^ 2 about zero energy shows the existence of 
a zero bias maximum (the denominator can be seen not to alter this property as long as the 
total transmission of the double barrier system is not too large, T v < 0.55). The zero bias 
anomaly shows up as a characteristic property of the overall conductance, see Figures 3 and 
4 (dotted lines). The zero bias maximum coincides with the maximum of the conductance 
product G (e) G (— e) of the corresponding NININ junction at zero energy. Figures 3 and 
4 illustrate the crossover from zero to finite bias anomalies for two different interbarrier 
distances d as the strength of barrier Ii is increased and I2 is kept fixed. For an interbarrier 
distance larger than the coherence length of the superconductor, several Andreev resonances 
show up (see Fig. 4). Note that the inversion of the barrier sequence transforms the zero 
voltage conductance from a local minimum to a local maximum, while keeping the same 
zero voltage conductance. This is illustrated in Fig. 4 by the pair of solid lines around zero 
voltage corresponding to inverse barrier sequences. 

In an attempt to understand the width of zero and finite bias peaks as well as the 
position of the first finite bias anomaly, we have compared them to the Thouless energy of 
the system. The Thouless energy!! E c in a disordered system can be defined as the product 
of the dimensionless normal state conductance g and the level spacing 5E, E c = gSE. We 
have determined a corresponding Thouless energy in our system by taking the level spacing 
of the INI box multiplied by its overall conductance. In the numerical simulations of a weakly 
transparent barrier system, we found good agreement of this energy scale with both the width 
and the position of the finite bias anomaly. The width of the zero bias anomaly coincided 
with the characteristic energy scale of the conductance correlator (G {E + e) G (E)) e as 
well as with E c . In this case, the distribution of the channel transmission of the double 
barrier system is bimodal and resembles strongly to that of a dirty systemll. As the overall 
transparency approaches 1, however, the width of the resonances and E c disagree. In this 
limit, the double barrier system does not represent the bimodal distribution of the dirty 
system properly and we cannot expect the above definition of the Thouless energy to be a 
sensible quantity. Of experimental interest is the effect of temperature on the conductance 
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resonances predicted above. A finite bias anomaly will be smeared out to a zero bias anomaly 
for a temperature of the order of the Thouless energy. In the systems investigated here this 
energy is of the order of 0.1A. 

The interference of multiple scattering processes between the scattering region and the 
NS interface thus produces an interesting structure in the differential conductance. Recently, 
experiments in disordered NS and SmS junctionJll^ have concentrated on the observation 
of zero and finite bias anomalies. It has been understood theoretically^HElil that these 
features are due to the interplay between the barrier at the NS interface and disorder in the 
normal lead. At small disorder, the differential conductance exhibits a zero bias maximum, 
while at large enough disorder, a finite bias peak is expectedB at a voltage of the order 
of the Thouless energy E c of the normal lead. This has recently been confirmed in an 
experimentEl. Here, we have found the existence of the analog features in a ballistic double 
barrier NS junction (in addition, the double barrier junction shows higher harmonics in 
the resonances of the conductance). The ballistic point of view applied to the disordered 
NS junction thus suggests the interpretation of the finite bias anomaly as a superposition of 
resonances due to quasi-bound Andreev states between the superconductor and the disorder. 



V. CONCLUSION 

We have used the scattering matrix approach including finite voltage and nonzero tem- 
perature dependences to describe the current- volt age characteristics of NS junctions. The 
current- volt age relation has been expressed through a spectral conductance which takes into 
account the explicit voltage dependence of the scattering problem. We have demonstrated 
the existence of a symmetry of the spectral conductance under sign reversal in the sub- 
gap regime. Above the gap, the symmetry is destroyed as the conductance approaches the 
Landauer expression. We have presented the spectral conductance formula in the Andreev 
approximation for a multichannel NS junction at all voltages in (|1"3|). This result has en- 
abled us to carry out a study of a double barrier NINIS junction at finite voltage, which has 
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revealed the crossover from normal to Andreev resonances in the INI interferometer as the 
ratio of the barrier strength is varied. Interestingly, a radical dependence of the conductance 
on the sequence of the barriers has been found. We have learned about the mechanism gov- 
erning the occurrence of and the crossover between zero and finite bias anomalies, tracing 
them back to the energy dependence of the modulus and the phase of the reflection am- 
plitude, respectively. Finally, we have established the similarity of the conductance of the 
double barrier NS junction with a disordered NS junction. 
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FIGURES 
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FIG. 1. Schematic structure of a disordered NS junction. Ballistic normal (iVi, N2) and su- 
perconducting leads (S) are coupled to reservoirs at chemical potential [i — eV and fx, respectively. 
Scattering is limited to the hatched region between the ballistic leads Ni and N2. 
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FIG. 2. Andreev resonances and resonance condition for the phase ip. Bottom: phase ip (e) 
of the reflection amplitude versus energy. The solid line represents the I1NI2 interferometer with 
barriers strengths Hi = 2%vp (R\ = 0.8) and H2 = Tivp (R2 = 0.5), the dashed line stands for the 
inverse barrier sequence (R\ = 0.5 and R2 = 0.8). The Andreev resonance condition for the phase 
is met for a pair of energies ±e n with phase difference Aip (e n ) = tt + 2nir. This phase condition 
can be fulfilled only by the first barrier sequence (R\ > R2, solid line) at those energies indicated 
in the graph, and produces a peak in the conductance of the NINIS structure. Interchanging the 
strengths of the two barriers, the Andreev resonances are greatly reduced, although the NININ 
double barrier system has the same overall transmission probability T (dotted line) as determined 
by the normal resonances. Top: conductance (arbitrary units) of the double barrier NS junction 
versus energy, the solid line again representing the barriers R\ = 0.8, R2 = 0.5, and the dashed 
line the barriers R\ = 0.5, R2 = 0.8. Note the symmetry of the resonances with respect to e = 0, 
which is due to their electron-hole character. 
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FIG. 3. Differential conductance for a multichannel NINIS junction with interbarrier distance 
d = 2vp/A = 2ir£. The average conductance per channel is plotted versus voltage at temperature 
T = 0. The barrier strengths H (reflection probabilities R = (1 + H 2 )/H 2 ) in units of Tivf are 
Hi = 0.5 (Ri = 0.2) (dotted line), 1.0 (0.5) (dashed), 1.5 (0.7) (dot-dashed), and 2.0 (0.8) (solid) 
in the normal region, while Hi = 1 (i?2 = 0.5) at the NS interface. The normal state conductance, 
which is roughly independent over this voltage range, is indicated on the left. With increasing 
barrier strength H\ the zero bias anomaly develops into a finite bias anomaly as the Andreev 
resonance is formed for H\ > H^. 
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FIG. 4. Differential conductance for a multichannel NINIS junction of width d = 4-uf/A = 47r£. 
The average conductance per channel is plotted versus voltage for temperature T = 0. The barrier 
strengths (reflection probabilities) in units of Tivf are H\ = 0.2 (i?i = 0.04) (dotted line), 0.5 
(0.2) (dashed), 0.8 (0.4) (dot-dashed), and 1.1 (0.54) (solid) in the normal region, and fixed at 
H2 = 0.5 (i?2 = 0.2) at the NS interface. The normal state conductivity is indicated on the left. 
As the strength of Hi is increased, the zero bias anomaly turn into an finite bias anomaly and 
several Andreev resonances appear. For the last choice of barriers Hi = 1.1 and H2 = 0.5 we 
have interchanged the barrier sequence: the conductance at zero voltage remains the same while 
changing from a zero bias minimum to a zero bias maximum (short solid line). 
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